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AERONAUTICAL SYMBOLS 
L FUNDAMENTAL AND DERIVED UNITS 



Length _ 
Time- 
Force- _ 



Symbol 



I 
t 

F 



Metric 



Unit 



meter 

second 

weight of one kilogram _ 



Svmbol 



m 
s 

kg 



English 



Unit 



foot (or mile) 

second (or hour) 

weight of one pound __ 

horsepower 

mi./nr 

ft. /sec 



Symbol 



ft. (or mi.) 
sec. (or hr.) 
lb. 



Power _ 
Speed. 



kg/m/s. 
fkm/h__ 
\m/s 



k. p. h. 
m. p. s. 



hp 

m. p. h. 
f. p. s. 



W, Weight == mg 

g, Standard acceleration of gravity = 9. S0665 
rn/s 2 = 32.1740 ft./sec. 2 
W 

m, Mass = — 
9 

p, Density (mass per unit volume). 
Standard density of dry air, 0.12497 (kg-m~ 4 b, 

s 2 ) at 15° C. and 760 mm = 0.002378 c, 

(lb.-ft" 4 sec. 2 ). 
Specific weight of "standard" air, 

kg/m 3 = 0.07651 lb./ft. 3 . 



2. GENERAL SYMBOLS, ETC. 

mk 2 , Moment of inertia (indicate axis of the 
radius of gyration k, by proper sub- 
script). 
Area. 

Wing area, etc. 
Gap. 
Span. 
Chord. 



S, 



1.2255 S' 



Aspect ratio. 
CoefFicient of viscosity. 



r, 



True air speed. 



Dynamic (or impact) pressure = ^pT 
Lift, absolute coefficient C L =^, 



3. AERODYNAMICAL SYMBOLS 

Q f Resultant moment. 
1 ' Resultant angular velocity. 



Drag, absolute coefficient C D = 



D 

qS 



D„, Profile drag, absolute coefficient Cj>^=-J 



D t , Induced drag, absolute coefficient C/> ( = 
D 



D 



Parasite drag, absolute coefficient C Dp —^ 
Ct Cross-wind force, absolute coefficient 

Cc ~qS 
R, Resultant force. 

i w , Angle of setting of wings (relative to 
thrust line). 

ii tJ Angle of stabilizer setting (relative to 
thrust line). 



yi 

p — > Reynolds Number, where I is a linear 

dimension. ■ 
e. g., for a model airfoil 3 in. chord, 100 
mi./hr. normal pressure, at 15° C, the 
corresponding number is 234,000; 
or for a model of 10 cm chord 40 m/s, 

the corresponding number is 274,000. 
Center of pressure coefficient (ratio of 
distance of c. p. from leading edge to 
chord length). 
a, Angle of attack, 
e, Angle of down wash. 
a 0 , Angle of attack, infinite aspect ratio. 
a iy Angle of attack, induced. 
a a , Angle of attack, absolute. 

(Measured from zero lift position.) 
7 Flight path angle. 



REPORT No. 413 



A METHOD FOR COMPUTING LEADING-EDGE LOADS 



By RICHARD V. RHODE and HENRY A. PEARSON 
Langley Memorial Aeronautical Laboratory 



93261—32 



1 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



NAVY BUILDING, WASHINGTON, D. C. 



(An independent Government establishment, created by act of Congress approved March 3, 1915, for the supervision and direction of the scientific 
study of the problems of flight, Its membership was increased to 15 by act approved March 2, 1929 (Public, No. 908, 70th Congress). It consists 
of members who are appointed by the President, all of whom serve as such without compensation.) 



Joseph S. Ames, Ph. D., Chairman, 

President, Johns Hopkins University, Baltimore, Md. 
David W. Taylor, D. Eng., Vice Chairman, 

Washington, D. C. 
Charles G. Abbot, Sc. D., 

Secretary, Smithsonian Institution, Washington, D. C. 
George K. Burgess, Sc. D., 

Director, Bureau of Standards, Washington, D. C. 
Arthur B. Cook, Captain, United States Navy, 

Assistant Chief, Bureau of Aeronautics, Navy Department, Washington, D. C. 
William F. Durand, Ph. D., 

Professor Emeritus of Mechanical Engineering, Stanford University, California. 
Benjamin - D. Foulois, Major (icneral, United States Army, 

Chief of Air Corps, War Department, Washington, D. C. 
Harry F. Guggenheim, M. A., 

The American Ambassador, Habana, Cuba. 
Charles A. Lindbergh, LL. D., 

New York City. 
William P. MacCracken, Jr., Ph. B., 

Washington, D. C. 
Charles F. Marvin, M. E., 

Chief, United States Weather Bureau, Washington, D. C. 
William A. Moffett, Rear Admiral, United States Navy, 

Chief, Bureau of Aeronautics, Navy Department, Washington, D. C. 
Henry C. Pratt, Brigadier General, United States Army, 

Chief, Materiel Division, Air Corps, Wright Field, Dayton, Ohic. 
Edward P. Warner, M. S., 

Editor "Aviation," New York City. 
Orville Wright, Sc. D., 

Dayton, Ohio. 



George W. Lewis, Director of Aeronautical Research. 
John F. Victory, Secretary. 
Henry J. E. Reid, Engineer in Charge, Langley Memorial Aeronautical Laboratory, Langley Field, Va. 
John J. Ide, Technical Assistant in Europe, Paris, France. 



EXECUTIVE COMMITTEE 



Joseph S. Ames, Chairman. 
David W. Taylor, Vice Chairman. 



Charles G. Abbot. 
George K. Burgess. 
Arthur B. Cook. 
Benjamin D. Foulois. 
Charles A. Lindbergh. 
William P. MacCracken, Jr. 



Charles F. Marvin. 
William A. Moffett. 
Henry C. Pratt. 



Edward P. Warner. 
Orville Wright. 



John F. Victory, Secretary. 



REPORT No. 413 



A METHOD FOR COMPUTING LEADING-EDGE LOADS 

By Richard V. Rhode and Henry A. Pearson 



SUMMARY 

In this report a jormula is developed that enables the 
determination of the proper design load for the portion of 
the wing forward of the front spar. The formula is inher- 
ently rational in concept, as it takes into account the 
most important variables that affect the leading-edge load, 
although theoretical rigor has been sacrificed for simplicity 



INTRODUCTION 

Recent failures of the leading-edge structures of 
some airplanes at high angles of attack and in nose 
dives have indicated the necessity for a revision of the 
specifications for the design of the leading edge of the 
wing. (See figs. 1 and 2.) While the present Army and 
Navy design rules (references 1 and 2) furnish a fairly 




Figure I. — Leading-edge failure that occurred in a fast dive 



and ease of application. Some empirical corrections, 
based on pressure distribution measurements on the 
"PW-9" and "M-8" airplanes, have been introduced to 
provide properly for biplanes. 

Results from the formula check experimental values in 
a variety of cases with good accuracy in the critical load- 
ing conditions . The use of the m et hod for design purposes 
is therefore felt to be justified and is recommended. 



good criterion for the strength of the leading-edge struc- 
ture in the high angle of attack condition, the entirely 
arbitrary provision for the nose-dive condition is 
inadequate in many cases. The Aeronautics Branch of 
t he Department of Commerce has no rule for either case. 

The National Advisory Committee for Aeronautics 
is now in a favorable position to study problems of this 
nature, making use of pressure-distrubution data from 

3 



4 REPORT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



fiighl testsand from tests at high Reynolds Numbers in 
the variable-density wind tunnel which have been accu- 
mulated over a period of years. 

It is, therefore, the purpose of this paper to develop 
and to present a more satisfactory rule for the practical 
determination of leading-edge design loads than has 
heretofore existed. 

The development of the formula involves, basically, 
theoretical considerations set forth in reference 3, 
although, in keeping with the interests of the practical 
designer, liberties have been taken with the rigorous 
theory to simplify the result as much as possible con- 
-1-1 ent with reasonable accuracy. Also, although the 
formula gives fair results for any condition of flight, 



DERIVATION OF THE FORMULA 

In reference 3, Theodorsen shows that the total load 
on a wing section may be considered as the sum of a 
"basic" load, which is a function only of the shape of 
the mean camber line, and an additional load, which is 
a function of the angle of attack measured with respect 
to the "ideal" angle at which the basic load occurs and, 
to a minor extent, a function of the nose curvature. 
He shows, further, that the distribution of the basic 
load is a function of the mean camber, and that the dis- 
tribution of the additional load is the same for any air- 
foil except for a narrow region near the leading edge 
where it becomes dependent upon the radius of the nose. 

Figure 3 illustrates these points. 




FlGTTBE 2— Result of leading-edge failure. Air entering at the opening of the seam burst the fabric 



il has been adjusted to give the best accuracy in the 
critical loading conditions. 

The principal result sought has been the value of t he- 
shear at the forward face of the front spar. Thus, the 
formula is developed to obtain this result and does not 
include provisions for the rational determination of the 
moment. An examination of pressure diagrams, how- 
ever, indicates that the centroid of the diagram area 
forward of the front spar, whatever the spar location, 
is confined within fairly narrow limits in terms of per- 
centage of spar distance from the leading edge for all 
airfoils. Empirical rules for the location of the center 
of gravity of the leading-edge load are, therefore, 
derived from which a static test may be devised to give 
a reasonably correct value of the moment as well as the 
shear at the critical section. 



Neglecting the minor variations at the nose, we may 
write, with only small error, for the whole wing section. 



C N = C B +K (a-a,) 



(1) 



where, 



Oi 
K 



-total load coefficient 
-basic load coefficient at a f 
-constant additional load per radian 
a — nominal or geometric angle of attack 

(radians) 
a 7 — ideal angle of attack 
For the present purpose it is of interest to examine the 
portion of the load forward of the front spar location, 
which is usually anywhere from 8 to 15 per cent of the 
chord. Equation (1) may be modified as follows: 



C s = k C B +K' (a- a,) 



(2) 
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where, C s is the load coefficient forward of front spar ar/° = 6.23 (2/1 — 2/5) +0.47 (2/2 — Z/0 

face and k and K' are appropriate constants. 

Equation (2) is in such a form that the leading-edge where y u y 2y 1/4, and >/ 5 are the ordinates of the mean 

load is given as a fraction of the total load on the sec- camber line, as fractions of the chord, with respect 

tion. The use of the function (a — ai) restricts the to a line joining the extremities of the mean camber 




> ct 1 



ct - or. 1 




V 



Ct < OLj 

Figure 3.—Pressure distribution above and below the "ideal" angle of attack 

application of the formula to two-dimensional flow, and, line, at stations ^ = 0.00542 c, r 2 =,0.125c, & 4 = 0.874 

for a practical case, the induced angle of attack a t c, and j 5 = 0.995 c. 

would have to be determined and subtracted from a. The formula may be simplified and its use facilitated 

Further, a 2 would have to be determined from the by the substitution of a function of C iV for («-«/). 

Gauss solution of Theodorscn's expression, viz: | Referring to Figure 4, 
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(a - a 7 ) « 



On Ob 
AC„ 
Aa 



. Ms rr 
where, — — = 5.5. 

Substituting in equation (2), the following expression 
is obtained: 



(3) 



At this point it is necessary to determine /r, C B , and 
K'. It has been pointed out that the distribution of 
the additional load which arises as a result of any depar- 
ture from a f is constant. The value of K', therefore, 
depends only upon the spar location, and graphical 
integration of the additional load-curve area forward of 



value of « — «/is also large. The argument does not 
apply to airfoils with very low camber near zero lift, 
but in such cases the total loads are small and of little 
interest. It is, therefore, assumed, for the sake of sim- 
plicity, that the value of C B is a function only of the 
general shape of the median line and of the maximum 
mean camber (measured always with respect to the 
chord of the mean camber line), and that the value of 
k y which is the portion of the basic load forward of the 
spar, is a function only of spar location for airfoils of the 
same general shape. To obtain working values of C B 
and k, curves have been drawn through theoretically 
derived values for several airfoils of conventional and 
reflex form with various cambers and, in the case of k, 
for several spar locations. 
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FlOUBE 4.— Relation between Cs and Cn 

any particular front spar face yields the desired result. 
C B is a function of the mean cambej line and varies 
I'roni zero for symmetrical sections to finite positive 
quantities for cambered airfoils. It may be considered 
negative for inverted sections. While its value may be 
determined precisely only by means of pressuie meas- 
urements at the ideal angle of attack or by Theodor- 
sen's expression, an approximation is sufficient for the 
present purpose. It may be said here, in justification 
of this step, that the ideal angle of attack is the angle 
at which the flow enters the leading edge smoothly, 
corresponding to the Kutta condition for smooth flow 
at the trailing edge. The basic load on the forward 
portion of the section is therefore small compared to 
the additional load imposed when the airfoil is at the 
angle of attack of maximum lift, one of the critical 
loading conditions for the leading edge. The same 
argument applies for the other critical condition, that 
of nose dive, in which, for commonly used airfoils, the 



















































Fo 


r- Oil O 


thers 




























































/ / 


^Fo 


~ refl 


exed 


lirfoi 


Is 



























































































































■ 



















O 5 10 15 20 

Location of forward face offronf spar in % of chord, fx) 

Figure 5. — K\ against spar location 

Before t lie formula is put in its final form, the normal 
discrepancy between theoretical and experimental re- 
sults must be taken into account. This discrepancy 
arises largely as a result of skin friction and is evidenced 
by progressively increasing pressure losses as the trail- 
ing edge is approached. The effect is, therefore, to 
shift the line of action of the experimental total load 
forward of that for the theoretical load, which results 
in an increase in the values of k and K' . The multi- 
plying factor for k and K f was found, by a .method of 
averages, to be 1.17. 

The formula may now be written 



a=i.i7 



1.17 K'C N 
5.5 



(4) 
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The result is further condensed to 

Cs=(-KiCb+K 2 C n ) (5) 

which is the final and useful form, giving what may be 
termed the " leading-edge shear coefficient, " or the 
Leading-edge load per unit chord, per unit span, per 
unit </. From this, substituting C L for C Ni the total 
load per unit span is 

w Le =CsY2pV 2 c 

= (-K 1 C B +K 2 C L )cxy 2P V 2 ) 

where c, p, and V have their usual significance. 

Curves for Ki and K 2 as functions of spar location 
are given in Figures 5 and 6, respectively, and curves of 

•^i 1 1 1 1 1 1 . 1 1 
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Location of forward face of front spar in 
per cent of chord, (x) 

FlGUBE 6.— K% against spar locution 

C B against maximum moan camber arc given in Figure 7. 
Figure 8 is included merely to show the manner in 
which y max is measured. 

COMPARISON OF RESULTS FROM FORMULA AND 
FROM PRESSURE-DISTRIBUTION TESTS 

For the purpose of checking the validity of the for- 
mula by comparison with pressure-distribution dia- 
grams, only those diagrams which were obtained from 
the variable-density wind tunnel at high Reynolds 
Numbers (reference 4) and from flight tests are used. 
While a vast quantity of pressure-distribution data 
from other sources is available, they have shown such 
inconsistency among themselves and with variable- 
density wind tunnel and flight results that it is 
believed advisable to avoid confusion, and possibly 



misleading conclusions, by the elimination of them 
from consideration altogether. 

A comparison of the calculated shear coefficients C s 
with experimental values obtained from tests on mono- 
plane airfoils in the variable-density tunnel is given in 
Table I. It will be noted that the agreement is good at 
high angles of attack for the variety of airfoils and spar 
locations given, the maximum difference being 12.9 per 
cent in the case of the N. A. C. A. 84-J airfoil with spar 
location at 20 per cent chord. 

At lift coefficients of zero or slightly below, represent- 
ing the nose-dive condition, the agreement is quite good 
at all spar locations for the R. A. F. 30 and N. A. C. A. 
84 airfoils. Larger errors, however, are apparent for 
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ymox i n P& r cent chord 

Figure 7.— Relation of Cn to maximum mean camber 




PlOUKE 8.— Maximum mean camber, j/,„» s of the Clark Y airfoil 



the M-6, Clark Y, and 84-J. There are several rea- 
sons for this apparent decrease in accuracy. First, 
differential normal pressures are less at the nose near 
zero lift than at maximum lift when the air speed is 
maintained constant. This condition results in larger 
percentage experimental errors at the low angles of 
attack, since small pressures are more difficult to meas- 
ure with accuracy than large pressures. Second, the 
peculiar shape of the pressure diagram for lift coeffi- 
cients near zero is of itself a cause for greater experi- 
mental error. Slight inaccuracies in locating pressure 
orifices result in large percentage errors in the leading- 
edge shear coefficient since the pressure gradient along 
the chord near normal front-spar locations is extremely 
steep. This source of error is an important considera- 
tion when using test data from small models. In 
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addition, since the leading-edge part of the pressure 
diagram is roughly a narrow triangle in shape (see fig. 
3), much depends upon the proper location of orifices 
and upon accurate measurement at a pressure station 
near the apex of the "triangle." Such causes as the 
above could easily explain most of the apparent 20 per 
cent error in the results on the Clark Y. 

In the case of the M-6 another source of error exists. 
Because this airfoil lias such a low camber, the value of 
a— a/is small near zero lift ; thus (he basic load, which 
in the formula is only approximated, largely predomi- 
nates. In fact, it may be expected that, as a — ai 
approaches such a value that the difference between 
the "basic" pressures and the "additional" pressures 
near the nose approach zero, percentage errors will 
become infinitely great. These errors, however, have 
no practical significance. It is to be noted, in this con- 
nection, thai the nose-dive shear coefficients for tin 4 
M 6 are quite small compared with bhose for the other 
airfoils, excepting the symmetrica] K. A. E\ 30 section. 

The greatest errors occur in the case of the NT. A. C. 
A. 84-J airfoil slightly below zero lift. Pressure dia- 
grams for this condition, showing definitely that the 
lower surface has stalled, indicate that these errors are 
almost entirely a result of the abnormal shape, or mean 
camber, of the section. It will be noted that the error 
is almost entirely eliminated if the concavity in the 
lower surface is removed to prevent this stalling; e. g., 
it will be noted that the error is small in the ease of the 
N. A. C. A. 84, which is simply the 84-J with the con- 
cave lower surface replaced by a flat lower surface. 

In view of the above discussion, the accuracy of the 
results obtained with the formula applied to mono- 
planes is considered reasonably good. It is to be 
doubted that better results could have been obtained 
without greater precision in test measurements and 
appreciably greater complications in the method or 
formula. 

In Table II results obtained with the formula are 
compared with experimental results obtained in flight 
on the M-3 and PW-9 airplanes having the Clark Y 
and Gottingen 436 sections, respectively. The test 
data represented in this table are of an appreciably 
higher order of accuracy than those in Table I. having 
been recently obtained after improvements in test 
methods were effected. The maximum error on the 
Clark Yin this table is -11 per cent, an amount which 
is not greater than might be expected as a result of the 
biplane arrangement. The upper wing in this case is 
alone represented . 

The error on the upper wing of the PW-9 airplane is 
about the same in magnitude and of the same sign as 
the error on the M-3 upper wing in the high angle of 
attack condition. In the nose-dive condition on the 
PW 9 lower wing, however, the error is consistently 
high, averaging about 30 per cent, while on the upper 
wing in this condition the error is negligible. 



EFFECT OF BIPLANE ARRANGEMENT 

No attempt is made here to include rationally tin 4 
effects of superposed wings as in biplane combinations. 
The character of these effects should, however, be 
understood and some provision made for them in the 
design. 

One airfoil mounted beneath another airfoil may be 
considered to have two effects on the latter. First, by 
virtue of its downwash at positive values of lift, tin 4 
lower wing causes the upper wing to operate at an 
effective angle of attack which is smaller than that at 
which it would operate as a monoplane. This 
effect does not influence appreciably the shape of the 
pressure diagram. Second, by virtue of the camber 
of the upper surface of the lower wing, the streamlines 
are curved even at an appreciable distance from the 
wing, a phenomenon which results in a decrease in the 
effective camber of the upper wing. This effeel causes 
a small forward shift of the center of pressure on the up- 
per wing at high angles of attack and an increase in the 
leading-edge load. Thus, it is to be expected that the 
formula for leading-edge loads will give low values for 
the upper wing at high angles of attack. On the M S 
and PW 9 airplanes (Table II) the result is 11 per cent 
too low. Most of this error, however, may be experi- 
mental error and error from the formula. In view of 
this and the reasonably small magnitude of the error, 
it is not considered necessary nor advisable to make any 
correction for the upper wing at high angles of attack. 

The effect of the upper wing on the lower is likewise 
small at the high angles of attack. 

At or near zero lift the effect of the curvat lire of i he 
streamlines appears to be small on the upper wings of 
both the M-3 and PW-9. On the PW 9 lower wing, 
however, the effect appears to be quite pronounced, as 
has been previously shown. A number of careful 
tests near zero lift on this airplane all bear out the fact 
that the shear coefficient on the lower wing is about 
30 per cent greater than that on the upper wing or that 
obtained from the formula. It is interesting to recall, 
in this connection, thai the leading-edge failure on a 
recently built diving bomber occurred on the lower 
wing in a dive. (Figs. 1 and 2.) This evidence, with 
the results of the PW-9 tests, indicates that the 
requirements for the lower wing of a biplane in the 
nose-dive condition should be increased over those 
for the upper wing or the monoplane. On the basis of 
the PW-9 tests, it is suggested that this increase should 
be 30 per cent in some cases. Since the effect is prob- 
ably caused by an induced change in camber resulting 
from curved air flow, it is more logical to include the 
effect in the formula by increasing the value of C B 
rather than by increasing the final result arbitrarily by 
30 per cent. From the PW-9 tests it appears that the 
value of C B should be increased 40 per cent for the 
biplane lower wing in the nose-dive condition. This 
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correction is recommended for general use until more 
information is available. 

USE OF THE FORMULA 

In any design, the strength of the leading edge should 
be investigated for the two critical conditions, (a) high 
angle of attack, and (b) nose dive. Experience, as well 
as analysis, has shown that no other condition need be 
considered and that in some cases, such as designs 
making use of symmetrical or nearly symmetrical air- 
foils and designs of large airplanes that are never dived, 
the high angle of attack condition only is of interest. 
To make use of the formula for these conditions, the 
leading-edge shear coefficient C s , and the dynamic 
pressure K pV 2 must be known. Cs, as has been shown, 
is a function of the shape of the airfoil and of the value 
of C N corresponding to the condition being inves- 
tigated. In a given design, with the airfoil and spar 
locations known, the constants K u K 2 , and C B are 
readily determined from the curves given in Figures 
5, 0, and 7. It is necessary to determine the values 
of C N and %pV 2 that will make the strength of the 
leading edge consistent with the strength of the rest of 
the wing structure. This may be done as follows: 

Case I — High angle of attack. — The primary wing 
structure is designed to fail in a condition correspond- 
ing to maximum C N with a certain specified load factor. 
From the general lift equation, this is equivalent to 
saving, 

nW=c Nm&x Sxy 2Po 17 

where, 7? IT — the load at failure 

n — high angle of attack load factor 
W — weight 
S — wing area 

p 0 — standard sea-level density 
Vi — the indicated speed at which, with Cv max , 
the equation is satisfied 
This speed Vi is the speed to use in the leading-edge 
formula for the high angle of attack condition, and 
6V max the proper value of C N . It is not essential that 
Cv max be determined with great precision, as any errors 
introduced in the shear coefficient will be approxi- 
mately compensated by errors of opposite sense in 
V 2 , with a resultant small error in the total leading- 
edge load. A representative problem, using the Clark 
Y airfoil, has indicated that C Nm&x may be in error by 
as much as 25 per cent and cause only a 5 per cent 
error in the total leading-edge load. 

In the case of biplanes, proper account should be 
taken of the relative wing-loading ratio. The above 
discussion of CAr max applies strictly only to monoplanes. 
It applies to biplanes when CV max is considered as that 
for the cellule and is used to determine T V. The mean 
lift coefficient for the biplane cellule should not be used 
in the formula to determine the leading-edge load for 
either upper or lower wings without correcting for the 
relative wing-loading ratio. This may be done by 



means of the following simple expressions, the wing- 
loading ratio being assumed known: 
L u + Lt = Wn 

or, 

L u = Wn-L x 
where, L u — lift on upper wing 
Li — lift on lower wing 
W — gross weight 
n — H. A. A. load factor 

Also, 

Tf_C N (upper) L u $ Wn -L K Si 
n C N (lowerpi, X S U ~ L t X S u 

Solve the above equation for L t . Then,, 
C N (lower) = y o 

and C N (upper) = RxC N (lower) where V? is the value 
found using the cellule C Nm&x and the high angle of 
attack load factor. 

As in the case of the monoplane, the biplane-wing 
lift coefficients found by the above method may not be 
true values, since they depend on V?, which itself has 
been found from an approximate cellule Cx m&x . How- 
ever, this makes no practical difference, any errors 
resulting in the shear coefficient being compensated by 
an error of opposite sense in V 2 to give a substantially 
correct leading-edge design load. 

The above biplane correction has nothing to do with 
the biplane corrections to the shear coefficient discussed 
in the preceding sections and is used merely to deter- 
mine the proper values of C N for the individual wings. 

Corrections to allow for the variation of C N along the 
span are not believed to be advisable in view of the 
added complication which would be involved. 

Case II — Nose dive. — In the nose-dive condition, 
the terminal velocity or the limited diving speed should 
be determined. The value of C N may be found by a 
solution of the conditions of static equilibrium for the 
case under consideration. It is suggested that, for the 
terminal-velocity dive, allowance be made for the possi- 
bility of encountering gusts and for slight inadvertent 
motions of the controls which may result in negative 
lift coefficients. This provision is important because 
the variation of leading-edge load with angle of attack 
near zero lift is extremely rapid, the load increasing 
greatly with small negative increments of lift coeffi- 
cient. So little is known about atmospheric conditions 
that it is difficult to establish a criterion for the de- 
termination of the proper negative lift coefficients on 
the basis of gusts. An examination of pilot-balloon 
data taken at Langley Field over a period of three 
months indicates that variations in horizontal wind 
velocities may be assumed as 15 feet per second, which 
would result in negative lift coefficients of from -0.15 
to —0.26 in the average case, depending on the termi- 
nal velocity. Other evidence exists which indicates 



10 REPORT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



that a value of 15 feet per second is not too conserva- 
tive. In addition to the possibilities of encountering 
gusts, there is also the fact that wings twist in dives 
under the heavy torsional moments experienced with 
the common wing sections. This influence may result 
in negative lift coefficients at the outer portion of the 
wing, even when the total lift coefficient is positive. 
Since this effect increases with speed, it is probably 
better to assume a constant negative lift coefficient for 
all cases instead of one which would vary approxi- 
mately inversely with the speed if a standard gust 
were used as a basis. Until more is known about 
conditions in the dive, it is felt that a value of C N 
not less than —0.2 should be used in the nose-dive 
analysis. 

The correction for relative wing-loading ratio is not 
to be used for the nose-dive analysis, but the 40 per 
cent increase in C n for the lower wings of biplanes, as 
recommended in the preceding section, should not be 
forgotten. 

TYPICAL PROBLEM 

( riven: 

Airplane Biplane pursuit 

Weight 2,720 pounds 

Area (upper) 184 square feet 

Area (lower) 88 square feet 

Area (total) 272 square feet 

Mean chord (upper) 5.75 feet 

Mean chord (lower) 4.00 feet 

Airfoil : Clark Y 

(2/max = 3.60 per cent; see fig. G.) 

II. A. A. load factor ft 12 

Spar-lace location x 10 per cent chord 

Relative wing-load ratio _ 1.2 

Terminal velocity 280 in. p. h. (410 

1'. p. s.) assumed 

Required : 

Shear at spar lace on both upper and lower 
wings in 

(a) H. A. A. condition 
(6) N. D. condition 

Solution: 

Constants: 
K x = 0.223 
K 2 = 0.367 
<7* = 0.525 
High angle of attack condition: 

Cvmax (cellule) = 1.4 (assumed) 
V2 _2nW 

2 X 12X2,720 „- 1AA ,, 
= 1.4X272X0.00237S = ' 2 ' in ° (f ' P " 

nW-L, & 

It — f * "o~ 

nW 12X2,720 . OAA , 
or Li — ~p q — r=7 t^t ^ = 9,300 pounds 

("1+0 ("*f +0 



9,300 



9,300 



1.23 



\ > 85.7X88 
and, CV M = 1.2X 1.23 = 1.48. 
Thus, for the upper wing: 

w l . e ={-K l C B +K 2 C N ) f %V?c u 

= ( - 0.223 X 0.525 +0.367 X 1 .48) 

X 72, 100X5.75 
= 210 pounds per foot span, or 



0.002378 



210 



; = 365 pounds per square foot average. 



0.002378 



0.1X5.75 

For the lower wing: 

Wi. e . = ( - 0.223 X 0.525 +0.367 X 1 .23) 

X72, 100X4 
= 118 pounds per foot span, or 

q i*x4 = 295 pounds per square foot average. 

The preceding values are the total design loads. 
Nose dive condition : 

Ki== 0.2231 

Ko = 0.367 as before 

C B = 0.525J 
Or.— -0.200 
C Nl = -0.200 



assumed for N. D. 



C B (corrected for lower wing) =1.4X0.525 = 0.735 
Thus, for the upper wing: 

fl 002378 

n.e. = ( - 0.223 X 0.525 - 0.367 X 0.200) 2 

X (410) 2 X5.75 
= 219 pounds per foot span. 

For the lower wing: 

0 00^ 7 S 

w Le . = (-0.223 X 0.735 - 0.367 X 0.200) "'V^ — 

X(410) 2 X4 
= 190 pounds per foot span. 

The above values are the total applied loads. The 
design load is obtained by multiplying by a factor of 
safety, say 2, which gives a result of 438 pounds per 
foot for the upper wing and 380 pounds per foot for the 
lower wing. 

It will be noted that it was not necessary to multiply 
the high angle-of-attack results by two, since the factor 
of safety was taken into account by using the design 
load factor in the determination of Vi 2 . The same 
result would have been obtained by calculating V? on 
the basis of the expected maximum applied load factor 
or n/2 (in this case, 6) and multiplying the final result 
by the factor of safety, 2. 

APPLICATION IN STATIC TESTS 

It is desirable, in static tests of the leading edge, to 
use a rectangular load distribution. Such a distribu- 
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tion is permissible if the stresses imposed at the critical 
section in flight can be well represented in this way. If 
the rectangular load distribution is to be used, the 
moment at the critical section (always the forward face 
of the spar) must, in addition to the shear, be approxi- 
mately correct. 

NO attempt has been made here to rationalize tin* 
determination of the moment of the leading-edge load 
about the face of the spar, since an empirical solution 
is believed to be within the limits of precision of prac- 
tical static tests. It has been found, from an examina- 
tion of a large number of pressure diagrams, that the 
location of the centroid of the part of the area forward 
of the front spar is, on the average, at 45 per cent of the 
spar location (0.45 x) in the high angle of attack condi- 
tion, and at 35 per cent (0.35 x) in the nose-dive 
condition. The relative position varies slightly with 
different spar locations and with different airfoils, but 
within the usual range the variation amounts to not 
more than three-fourths of 1 per cent of the total chord 
or only a small fraction of an inch lor the ordinary 
airplane. 

Static tests may, therefore, be made using a rec- 
tangular load distribution, the center of gravity of the 
load being at, 

0.45 x for II . A. A. 
0.35 x for N. D. 



Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., January 10, 1031. 
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TABLE I 

comparison of shear coefficients from 
formula and from pressure-distribution 
tests in the variable-density tunnel 

r. a. f.-30 airfoil 



Maximum mean camber=0% c; Cb=0 



Spar 
loca- 
tion 

(%c) 




Ki 


K 2 


Cs 


Error, 

% 




Cl 


(from 
curve) 


(from 
curve) 


Com- 
puted 


Experi- 
mental 


Remarks 


5 


±1.20 
±. 95 
±. 76 
±. 57 
±. 34 


0. 162 
. 162 
. 162 
. 162 
. 162 


0. 222 
.222 
.222 
. 222 
.222 


0. 26(5 
.211 
. 169 
. 127 
. 0755 


0. 253 
.224 
. 176 
. 117 
.074 


5. 1 
-5.8 
-4.0 
8. 5 
2.0 


II . A. A. condition. 


I () 


±1.20 
±. 05 
±. 70 
±. 57 
±.34 


.223 
.223 
.223 
. 223 
.223 


.361 
.361 
.361 
.361 
. 361 


.433 
.343 
.274 
.206 
. 123 


.443 
.376 
.294 
. 204 
. 129 


-2.3 
-8.8 
-6.8 

LO 

-4.7 


II. A. A. condition. 


15 


±1.20 
±. 95 
±. 76 
±.57 
±.34 


.252 
.252 
. 252 
.252 
.252 


.477 
.477 
.477 
.477 
.477 


. 572 
. 453 
.363 
.272 
. 162 


.589 
. 4X0 
.385 
.264 
. 171 


-2.9 
-6.8 
-5.7 
3.0 
-5.3 


II. \ . A . condition. 


20 


±1.20 
±. 95 
±. 76 
±.57 
±. 34 


.236 
. 236 
. 23(5 
. 236 
.236 


. 571 
.571 
.571 
.571 
.571 


.685 
..543 
.434 
. 326 
. 194 


.702 
.577 
. 453 
. 315 
.208 


-2.4 
-5.9 
-4.2 
3.5 
-6.7 


H. A. A. condition. 






N. A. C. A. M-6 AIRFOIL 






Maximum mean camber = 2.215% c; Cb— 


0.226 


Spar 
loca- 
tion 
W c) 




Ki 


A' 2 


Cs 


Error, 

% 




Cl 


(from 
curve) 


(from 
curve) 


Com- 
puted 


Ex peri- 
mental 


Remarks 




5 


1.08 
1.05 
-. 10 


0. 149 
. 149 
. 149 


0. 222 
.222 
.222 


0.206 
. 199 
-.0559 


0. 191 
.203 
-. 0639 


7.9 
-2.0 
-12.5 


jfl. A. A. condition. 

Approximately N. 
D. condition. 


10 


1.08 
1.05 
-. 10 


.192 
. 192 
. 192 


.361 
.361 
.361 


.347 
. 888 
-. 0795 


.337 
.345 
-.0808 


3.0 
-2.6 
-1.6 


Jll. A. A. condition. 

Approximately N. 
1). condition. 


15 


1.08 
1.05 
-. 10 


.192 
. 192 
. 192 


.477 
.477 
.477 


.472 
.458 
-.0911 


.448 
.459 

-.(K« 


5.4 
-.2 
9.4 


Jh. A. A. condition. 

Approximately N. , 
D. condition. 


20 


1.08 
1.05 
-.10 


167 
. 167 
.167 


.571 
.571 
.571 


.579 
.562 
-.0948 


.541 

.559 
-. 0782 


7.0 
.5 
21.2 


jll. A. A. condition. 

Approximately N. 
D. condition. 
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TABLE I— Continued 

COMPARISON OF SHEAR COEFFICIENTS FROM 
FORMULA AND FROM PRESSURE-DISTRIBUTION 
TESTS IN THE VARIABLE-DENSITY TUNNEL— Con. 

CLARK Y AIRFOIL 



Maximum mean camber=3.6%c; Cb = 0.525 



Spar 
loca- 
tion 

(%c) 


Cl 


A'i 
(from 
curve) 


Ki 
(from 
curve) 


Com- 
puted 


Os 

Experi- 
mental 


Error, 

% 


Remarks 


5 


1.49 
1.44 
.02 

-. 19 


0. 1G2 
. 162 
. 162 

. 162 


0. 222 
.222 
.222 

.222 


0.246 
.235 
-. 0805 

-. 127 


0.247 
.209 
-. 0799 

-. 103 


-0.4 
12.4 
.8 

23.3 


Jll. A. A. condition. 

Approximate zero 
lift. 

ApproximateN. D. 
condition. 


10 


1.49 
1.44 
.02 

-. 19 


.223 
.223 
.223 

. 223 


.361 
.361 
.361 

.361 


.421 
. 403 
-. 110 

-. 186 


.434 
.382 
-. 120 

-. 153 


-3.0 
-5.5 
-8.3 

21.6 


Jll. A. A. condition. 

Approximate zero 
lift. 

ApproximateN. 1). 
condition. 


15 


1.49 
1.44 
.02 

-. 19 


.252 
. 252 
.252 

. 252 


.477 
.477 
.477 

.477 


. 579 
. 556 
-. 123 

-.223 


.589 
. 529 
-. 137 

-. 190 


-1.7 
4. 9 
-10.2 

17.4 


J II. A. A. condition. 

Approximate zero 
lift. 

Approximate N. 1). 
condition. 


20 


1.49 
1.44 
.02 

-. 19 


. 236 
. 236 
.236 

. 236 


.571 
.571 
. 571 

. 871 


. 727 
.698 
-. 112 

-.232 


. 723 
. 657 
-. 140 

-.224 


.0 
6.2 
-20.0 

3.0 


Jll. A. A. condition. 

Approximate zero 
lift. 

Approximate N. 1). 
condition. 


N. A. C. A. 84 AIRFOIL 

Maximum mean camber =5.25% c; Cb=0.7G2 


Spar 
loca- 

t ion 

(%e) 


Cl 


Ai 

(from 
curve) 


A'2 

(from 
curve) 


Com- 
puted 


Cs 

Experi- 
mental 


E rror. 
% 


Remarks 




1.51 
1.37 
0 

-. 19 


0. 1(12 
. 162 
. 162 
. 162 


0. 222 
.222 
.222 
.222 


0. 212 
.181 
-. 123 
-. 166 


0.216 
. 180 
-. 120 
-. 159 


-1.9 
.0 
2.5 
4.4 


Jll. A. A. condition. 
Zero lift. 

Approximate N. 1). 
condition. 


10 
15 


1. 51 
1.37 
0 

-. 19 


.223 
.223 
.223 
.223 


.361 
.361 
.361 
.361 


.375 
.325 
-. 170 
-.239 

.528 
.462 
-. 192 
-.283 


.396 
.336 
-. 181 
-.239 


-5.3 
-3.3 
-6. 1 
0 


J II. A. A. condition. 
Zero lift. 

Approximate N. D. 
condition. 


1.51 
1.37 
0 

-. 19 


.252 
.252 
. 252 
.252 


.477 
.477 
.477 
.477 


.539 
.466 
-. 200 
-.284 


-2.0 
-.9 

-6.8 
-.4 


Jll. A. A. condition. 
Zero lift. 

Approximate N. D. 
condition. 


20 


1.51 
1.37 
0 

-. 19 


.236 
.236 
.236 
.236 


.571 
.571 
.571 
.571 


.682 
.602 
-. 180 
-. 289 


.669 
.585 
-.214 
-310 


1.9 
2.9 
-15.9 
-6.« 


Jll. A. A. condition. 
Zero lift. 

Approximate N. D. 
condition. 
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TABLE I— ContiniKMl 

COMPARISON OF SHEAR COEFFICIENTS FROM 
FORMULA AND FROM PRESSURE-DISTRIBUTIOM 
TESTS IN THE YAH I ABLE-DENSITY TUNNEL— Con. 

N. A. C. A. 84-J AIRFOIL 



Maximum mean camber=7.3% c; Cb=0.97 



Spar 
loca- 
tion 
(%c) 


Cl 


A'i 
(from 
curve) 


A' 2 
(from 
curve) 


Cs 


Error, 

% 


Remarks 


Com- 
puted 


Experi- 
mental 


5 


1.72 
1.63 
0 

-.16 


0. 162 
. 162 
. 162 
.162 


0. 222 
.222 
.222 
.222 


0. 225 
.205 
-. 157 
-. 193 


0.230 
.188 
-.163 
-.121 


—2. 2 
9.0 
-3.7 
59.5 


Jll. A. A. condition. 
Zero lift. 

Approximate N. D. 
condition. 


10 


1.72 
1.63 
0 

-.16 


.223 
.223 
.223 
.223 


.361 
.361 
.361 
.361 


.405 
.372 
-.216 
-.274 


.410 
.350 
-.242 
-.206 


-1.2 
6.3 
-10.7 
33.0 


Jll. A. A. condition. 
Zero lift. 

Approximate N. D. 
condition. 


15 


1.72 
1.63 
0 

-. 16 


.252 
.252 
.252 
.252 


.477 
.477 
.477 
.477 


.576 
.533 
-.244 
-.321 


.564 

. INS 

-. 286 
-.256 


2*1 
9.2 
-14.7 
25.4 


J II. A. A. condition. 
Zero lift. 

Approximate N. D. 
condition. 


20 


1.72 
1.63 
0 

-. 10 


.236 
. 236 
. 236 
.236 


.571 
.571 
.571 
.571 


.753 
.702 
-.229 
-.320 


.713 
. 022 
-.310 
-.288 


5.6 
12.9 
-26. 1 
11. 1 


Jll. A. A. condition. 
Zero lift. 

Approximate N. D. 
condition. 



TABLE II 

COMPARISON OF LEADING-EDGE LOADS AS COM- 
PUTED. TO THOSE LOADS OBTAINED IN ACTUAL 
FLIGHT 

CLARK Y AIRFOIL (UPPER WING DOUGLAS M-3) 



Maximum mean camber 3.6% c; Cb=0.525; Chord 5.667 feet 



Spar 




AT, 
vfrom 
curve) 


Kt 
(from 
curve) 










loca- 
tion 

(%c) 


Cnp 


Speed 
f. p. s. 


(Com- 
puted) 


(Exper- 
imen- 
tal) 


Error. 
% 


Remarks 


7.2 


1.533 

.0144 
-. 2344 


0. 191 

. 191 
. 191 


0. 285 

.285 
.285 


107.6 

128.0 
122.5 


26.3 

10.0 
10.9 


29. G 

9.0 
10.0 


-11.0 

10. 5 
5.6 


H. A. A. condi- 
tion. 



GOTTINGEN 436 AIRFOIL (UPPER WING BOEING 
PW-9) 

Maximum mean camber 3.55% c; Ca=0.515; Chord 5. 42 feet 





1.743 


0. 252 


0. 477 


144 


93.7 


103. 6 


- 9. 0 






1.649 


.252 


.477 


214 


193.7 


202.0 


-4.1 






1.682 


.252 


.477 


183 


145. 1 


161.7 


-10.3 


Ih. A. A. condi- 




1.810 


.252 


.477 


180 


153. 1 


165.3 


-7.4 


tion. 




1.660 


.252 


.477 


212 


191.7 


203.0 


-5.6 




15 


.0564 


.252 


.477 


360 


-86.0 


-73.2 


17.5 


Approximate 
zero lift. 




-. 3466 


.252 


.477 


211 


-84.6 


-85.6 


-1.2 






-.4313 


.252 


.477 


207 


-92.8 


-92.8 


0 






-. 3668 


.252 


.477 


205 


-82.6 


-83.6 


1.2 






-. 3795 


.252 


.477 


212 


-90. 1 


-89. 6 


.6 






.0401 


.252 


.477 


381 


-103.8 


-107.0 


-3.0 





( ; OTTINGEN 436 AIRFOIL (LOWER WING BOEING 

PW-9) 

Maximum mean camber 3.55% c; C B =0.515; Chord 4.23 feet 





1.298 


0. 2405 


0.415 


144 


43. 


3 


45.6 


-5.0 






1. 194 


.2405 


.415 


214 


85. 


7 


87.2 


-1.7 






1.257 


.2405 


.415 


183 


67. 


1 


70.4 


-4.7 




1.324 


. 2405 


.415 


180 


69. 


3 


72. 8 


-4.8 


tion. 




1.207 


.2405 


.415 


212 


85. 


2 


85.4 


-.2 






.0403 


.2405 


.415 


360 


-69. 


8 


-117.0 


-40.3 


Approximate 
zero lift. 




-.314 


.2405 


.415 


211 


-56. 


9 


-76.0 


-25.1 




12.2 


-. 358 


.2405 


.415 


207 


—58. 


9 


-76.8 


-23.3 






-.291 


.2405 


.415 


205 


-51. 


8 


-63. 2 


-18.0 






-.313 


.2405 


.415 


212 


-57. 


4 


-72.0 


-20.3 






.0378 


.2405 


.415 


366 


-72. 


8 


-100.0 


-27. 2 






.087 


.2405 


.415 


367 


-59. 


7 


-84.0 


-28.9 






. 0997 


.2405 


.415 


370 


-57. 


2 


-81.2 


-29.6 






.0913 


.2405 


.415 


363 


-57. 


0 


-78.0 


-26. 9 






. 130 


.2405 


.415 


361 


-45. 


9 


-66. 5 


-31.0 
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Positive directions of axes and angles (forces and moments) are shown by arrows 



Axis 


Force 
(parallel 
to axis) 
symbol 


Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


Designation 


Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal 

Lateral 


X 
Y 
Z 


X 
Y 
Z 


rolling 


L 
M 
N 


Y i. Z 

Z > X 

X > Y 


roll 


e 


u 

V 

w 


V 
Q 


pitching 

\ awing 


pitch 


Normal 


yaw 









Absolute coefficients of moment 
L 



qbS 



_ M 
Vm qcS 



0,= 



Angle of set of control surface (relative to neu- 
N tral position), 6. (Indicate surface by proper 

subscript.) 



4. PROPELLER SYMBOLS 



P, 



Diameter. 
Geometric pitch. 
p/Dj Pitch ratio. 
V, Inflow velocity. 
V S) Slipstream velocity. 



T, Thrust, absolute coefficient C T = — jtm 
Q, Torque, absolute coefficient ^Q = ~^ijJ^ 



Pj Power, absolute coefficient C P = 

Cs, Speed power coefficient • = -y -p— 

77, Efficiency. 

n, Revolutions per second, r. p. s. 

3>, Effective helix angle ~ tan" 1 (j>~;~) 



5. NUMERICAL RELATIONS 



I' hp = 76.04 kg/m/s = 550 lb./ft./sec. 
1 kg/m/s = 0.01315 hp 
1 mi./hr.= 0.44704 m/s 
1 m/s = 2.23693 mi./hr. 



1 lb. = 0.4535924277 kg. 

1 kg = 2.2046224 lb. 

1 mi. = 1009.35 m = 5280 ft. 

1 m = 3.2808333 ft. 



